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We discuss different resummations of large logarithms that arise in hard-scattering cross 
sections of quarks and gluons in regions of large and small x. The large-x logarithms are 
typically dominant near threshold for the production of a specified final state. These soft 
and collinear gluon corrections produce large enhancements of the cross section for many 
processes, notably top quark and Higgs production, and typically the higher-order correc- 
tions reduce the factorization and renormalization scale dependence of the cross section. 
The small-x logarithms are dominant in the regime where the momentum transfer of the 
hard sub-process is much smaller than the total collision energy. These logarithms are im- 
portant to describe multijet final states in deep inelastic scattering and hadron colliders, and 
in the study of parton distribution functions. The resummations at small and large x are 
linked by the eikonal approximation and are dominated by soft gluon anomalous dimen- 
sions. We will review their role in both contexts and provide some exphcit calculations at 
one and two loops. 



1 Introduction 

Particle physics in high-energy hadron colliders depends crucially on our ability to calcu- 
late cross sections to an ever increasing theoretical accuracy, which is achieved by the in- 
corporation of higher-order coiTcctions. Hard-scattering cross sections in perturbative QCD 
obey factorization theorems [1] that play a key role in the calculation of these corrections. 
Typically, the cross section for a process involving the collision of two hadrons (proton- 
antiproton at the Fermilab Tevatron or proton-proton at the CERN LHC) into a specified 
final state can be described as a convolution of non-perturbative parton distribution func- 
tions that describe the parton content of the hadron, and a partonic cross section that can be 
calculated order-by-order in perturbation theory. The short-distance partonic cross section 
involves the scattering of quarks and gluons. The partonic processes are of the form 

fliPl) + f2{p2)^F{p) + X, (1) 
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where /i and /2 represent partons (quarks or gluons), F represents an observed system in 
the final state, such as a top quark or a jet or a Higgs boson, and X represents any additional 
final-state particles. The factorization is described schematically by 

where ah-^h2^F is the physical cross section (total or differential) for the production of final 
state F in the scattering of hadrons hi and /12, (^fi/h^ is the distribution function for parton 
fi with momentum fraction xi of hadron hi, and (7/i/2^f is the partonic cross section. 
The coUinear singularities are factorized in a process-independent manner and absorbed 
into the parton distribution functions which are dependent on the factorization scale jip- 
The physical cross section is in principle independent of the factorization scale ji p and the 
renormalization scale /i/j, but in practice there is a strong dependence because we truncate 
the infinite perturbative series at finite order (typically next-to-leading-order (NLO) or next- 
to-next-to-leading-order (NNLO) in the strong couphng ag). The parton-level cross section 
expUcitly involves the standard kinematical invariants, s = (pi + P2)^> t = {Pi — 
u = {p2 — p)^, formed from the 4-momenta of the particles in the hard scattering. 

Near threshold, i.e. when the energy of the incoming partons is just sufficient to produce 
a final state without additional radiation, the production cross section receives significant 
corrections from large-x logarithms [2-5]. These logarithms arise from incomplete can- 
cellations between virtual terms and terms that describe soft-gluon emission. Since near 
threshold any additional radiation has to be soft, the large-x logarithms are especially im- 
portant in that kinematical region. Large-x resummation depends critically on the color 
structure of the process [4,6-8] as well as the kinematics [4,9]. 

Small-x logarithms arise when the perturbative scales characterizing the hard- 
subprocess are much smaller than the total collision energy. In this case resummation of 
logarithms of the form ln(l/x) becomes important. When the transverse scales of the out- 
going scattered particles are similar and large this resummation can be described by the 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation [10-14]. This equation is a lin- 
ear integral equation which leads to an exponential rise of the cross section. The slope 
of this rise can be interpreted as a perturbative construction of the QCD Pomeron. This 
Pomeron is considered the mediator of many QCD diffractive processes, such as diffractive 
vector meson production. 

The conditions by which the BFKL evolution should be valid are satisfied by jet produc- 
tion with large rapidity gaps. Phenomenological studies of this process with the summation 
of the terms ctg ln"(l/x) (leading-order kernel) are not very predictive since the value of 
the coupling is a free parameter and the Regge energy scale, a sort of factorization scale 
at high energies, can only be fixed at higher orders. Inclusion of the next-to-leading order 
corrections, a^^^ ln"(l/x) [15, 16], brings the predictions in closer agreement with data. 

In the next section we discuss large-x resummation and finite-order expansions of the 
resummed cross section through next-to-next-to-next-to-leading order (NNNLO). In Sec- 
tion 3 we present some applications of large-x resummation to various hard-scattering pro- 
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cesses, namely top-antitop pair production, single top quark production, T^-boson produc- 
tion at large transverse momentum, and Higgs boson production via hb H. In Section 4 
we present typical one-loop and two-loop calculations in the eikonal approximation that are 
needed in resummations for processes with massive quarks, such as heavy quark pair pro- 
duction. Section 5 discusses small-x resummation and applications of BFKL. We conclude 
in Section 6. 

2 Large-x resummations 

Large-x resummations depend crucially on the kinematics and color structure of the process 
under study. In single-particle-inclusive (IPI) kinematics we identify one particle F with 
momentum p. In pair-invariant-mass (PIM) kinematics we identify a pair of particles (such 
as a heavy quark-antiquark pair) with invariant mass squared Q'^. 

In general, the partonic cross section a includes soft corrections in the form of plus 
distributions Vi {xth) with respect to a kinematical variable that measures distance from 
threshold, with / < 2n — 1 at n-th order in ag beyond the leading order. In IPI kinematics, 
Xth is usually denoted as S4 (or S2) and is defined by S4 = s + t + u — "^^w?, where the 
sum is over the squared masses of all particles in the process. At threshold, S4 = 0. The 
plus distributions are then of the form 



where M is a hard scale relevant to the process, for example the mass m of a heavy quark or 
the transverse momentum of a jet. The plus distributions are defined through their integral 
with the parton distribution functions by 



In PIM kinematics, xt^ is usually denoted as 1 — x or 1 — z, with z = Q'^/s — > 1 at 
threshold. Then the plus distributions are of the form 



W(s4/M2) 



(3) 
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(6) 



4 



Nikolaos Kidonakis, Agustm Sabio Vera, Philip Stephens 



The highest powers of these distributions in the nth-order corrections are the leading loga- 
rithms (LL) with I = 2n ~ 1, the second highest are the next-to-leading logarithms (NLL) 
with I = 2n — 2, etc. (note that the counting of logarithms is different in the exponent and 
in the fixed-order expansions). These logarithms can be resummed in principle to all orders 
in perturbation theory. 



2.1 Exponentiation 

The resummation of threshold logarithms is performed in moment space. By taking mo- 
ments, divergent distributions in 1 — z (or 54) produce powers of In N, with N the moment 
variable: 



Jo [l — zj. m + l V ^ 



(7) 



If we define moments of the partonic cross section by d-{N) = J dz z^ ^<5"(^) (PIM) or 
by a{N) = J {dsA/s) e-^''''/^a{sA) (IPI), then the logarithms of N that appear in d{N) 
exponentiate. 

The resummation follows from the factorization properties of the cross section. We 
begin the derivation of the resummed cross section by first writing a factorized form for the 
moment-space infrared-regularized parton-parton scattering cross section, aj^f^^pi^,^), 
which factorizes as the hadronic cross section 

(^fif2-*F{N,e) = 4>f^/f^{N,iiF,e) 4)f^/f^{N, iiF,e) af^f^^piN, fip, m) , (8) 

with the moments of (p given by <^(A^) = Jq dx x^~^(p{x). We factorize the initial-state 
collinear divergences, regularized by e, into the parton distribution functions, 0, which are 
expanded to the same order in ag as the partonic cross section, and we thus obtain the 
perturbative expansion for the infrared-safe partonic short-distance function a. 

The partonic short-distance function a still has sensitivity to soft-gluon dynamics 
through its N dependence. We then refactorize the moments of the cross section as [4, 8] 

<^fif2^F{N,e) = i>f^/f^ {N,HF,€) i)f^/f^ {N,iJ,F,e) 

X hII^^-^ {asii^n)) Sif-^ W Jj {N, mf, e) + 0(l/iV) , (9) 

where are center-of-mass distributions that absorb the universal collinear singularities 

from the incoming partons, ifyx are A^-independent hard components which describe the 
hard-scattering, Sli is a soft gluon function associated with non-collinear soft gluons, and 
J are functions that absorb the collinear singularities from massless partons, if any, in the 
final state. 

H and S are matrices in color space and we sum over the color indices / and L that 
describe the color structure of the hard scattering. The hard-scattering function involves 
contributions from the amplitude of the process and the complex conjugate of the ampUtude, 
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HjL = hi. The soft function Sli represents the coupling of soft gluons to the partons 
in the scattering. The color tensors of the hard scattering connect together the eikonal lines 
to which soft gluons couple. One can construct an eikonal operator describing soft-gluon 
emission and write a dimensionless eikonal cross section, which describes the emission of 
soft gluons by the eikonal lines [4, 6-8]. 

Comparing Eqs. ([Hi and Q, we see that the moments of the short-distance partonic 
cross section are given by 

(Tf,f,^F{N,HF,f^R) = -J — -J — -H\j^ (asM) 

X S{t^^ (^^^a,{^^R)^X{JJ {N,^iF,e) . (10) 

All the factors in Eq. (ITOl ) are gauge and factorization scale dependent. The constraint that 
the product of these factors must be independent of the gauge and factorization scale results 
in the exponentiation of logarithms of N imp/4> and Sli [4,5]. 

The soft matrix Sli depends on N through the ratio M/ (Nfip), and it requires renor- 
malization as a composite operator. Its A^-dependence can thus be resummed by renor- 
malization group analysis [17-20]. However, the product HilSli needs no overall renor- 
malization, because the UV divergences of Sli are balanced by those of Hil- Thus, we 
have [4, 8] 



H'lL = n^r^ {^s') ,^Hcn\{zl)-' 



i=a,b 

4) 



DL 



Sli = {Zs)lbSbaZs,ai, (H) 



where and 5*^ denote the unrenormalized quantities, Zi is the renormalization constant 
of the ith incoming partonic field, and Z5 is a matrix of renormalization constants, which 
describe the renormalization of the soft function. Z5 is defined to include the wave function 
renormalization necessary for the outgoing eikonal lines that represent any heavy quarks. 

From Eq. (ITT]) , we see that the soft function 5^/ satisfies the renormalization group 
equation [4, 6-8] 

d d \ 

fJ-TT- + Pids)-^ ) Sli = -(.Ts)lbSbi - Sla{Ts)ai , (12) 



dfi dg, 



where /3 is the QCD beta function and = Airag. Ts is an anomalous dimension ma- 
trix that is calculated in the eikonal approximation by expUcit renormalization of the soft 
function. In a minimal subtraction renormalization scheme and with e = 4 — n, where n is 
the number of space-time dimensions, the soft anomalous dimension matrix is given at one 
loop by 



s (^s) = - — ^ReSe^o^s(5s,e) . (13) 
2 do.. 
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The process-dependent matrices Ts have been calculated at one loop for all 2 ^ 2 partonic 
processes; a compilation of results is given in [8]. In processes with trivial or simple color 
structure Ts is simply a function (1x1 matrix) while in processes with complex color 
structure it is a non-trivial matrix in color exchange. For quark-(anti)quark scattering, Ts 
is a 2 X 2 matrix [4, 18]; for quark-gluon scattering it is a 3 x 3 matrix [7]; for gluon- 
gluon scattering it is an 8 x 8 matrix [7]. Complete two-loop calculations of soft anomalous 
dimensions for processes with massless quarks have appeared in [21]. Selected two-loop 
results for heavy quark production appeared in [22]. We present a sample one-loop calcu- 
lation in Section 4.1 and a sample two-loop calculation in Section 4.2, both with outgoing 
massive quarks (see [23]). 

The exponentiation of logarithms of N in the ratios i/^/cj) and in the functions J in 
Eq. (fTOl) . together with the solution of the renormalization group equation (fT2l) . provide 
us with the complete expression for the resummed partonic cross section in moment space 
[4,6-8,24,25] 



a'-^'iN) = exp 



exp 



X exp 



' ^ -ff^/f^{as{n)) exp 2dc, f ^/3(a,(/i)) 



^T^I^Hf^f^-*^ {asim)) exp 
X 5-/'i-^2~*^ I Q,_, I ^ I I exp 




2^F 



(as(^)) 



(14) 



The sums over i = 1,2 run over incoming partons. The sum over j is over massless partons, 
if any, in the final state at lowest order. The resummed expression is valid for either IPI 
or PIM kinematics. In IPI kinematics Ni = N{—ti/M'^), where tj denotes t or u, and 



Nj = N{s/M'^), while in PIM kinematics Ni 



Nj = N. Also N = Ne"^, with 7^ the 



Euler constant. 

The first exponent in Eq. ([14)) arises from the exponentiation of logarithms of N in the 
ratios V'/'/' of Eq. (ITOl) . and is given in the MS scheme by 



dz 



1 - z 



dX 

/(l-^)2 A 



(15) 



with Ai{as) = En=Ms/TTTAf'\ At one loop, A. 



Ci which is C'f = {N^ 



1) / {2Nc) for a quark or antiquark and Ca = Nc for a gluon, with Nc the number of colors, 
while Af^ = CiK/2 with K = (7^(67/18 - C2) - 5n//9 [26], where rz/ is the number of 



quark flavors and C2 = ttVG. Also Vi = E^=i(as/7i")''z^f"'', with u\"> = d. 

The second exponent in Eq. (fT4l) arises from the exponentiation of logarithms of in 



,(1) 
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the functions Jj of Eq. (iTOl ). and is given by 

E'f^iNj) = dz— / A,{as{\s))-Bj[as{{l-z)s) 

Jo V — z I J(i-2)2 A 



Q,((l-z)2s)]} . (16) 



Here Bj = T.n=i{(^s/T^TBf' with B) ' equal to 3Cf/4 for quarks and /?o/4 for gluons, 
where (Sq = (IICa — 2nj)/3 is the lowest-order (i function. 

The third exponent in Eq. (fT4l ) controls the factorization scale dependence of the cross 
section, and 7/,//. is the moment-space anomalous dimension of the MS density 4>f,/fr 
The f3 function in the fourth exponent controls the renormalization scale dependence of the 
cross section. The constant takes the value k if the Bom cross section is of order a^. 
Explicit expressions for the functions in these four exponents, and related references, are 
assembled for convenience in Appendix A of Ref. [27]. 

As noted before, both H and S are process-dependent matrices in color space and thus 
the trace is taken in Eq. ([141 ). At lowest order, the trace of the product of H and S reproduces 
the Born cross section. The evolution of the soft function 5 follows from its renormalization 
group equation, ([T2l ). and is given in terms of the soft anomalous dimension matrix Ts- 

2.2 NNNLO expansions 

The exponentials in the resummed partonic cross section can be expanded to any fixed order 
in as and then inverted to momentum space to provide explicit results for the higher-order 
corrections. A fixed-order expansion avoids using a prescription to regulate the infrared 
singularities in the exponents and thus no prescription is needed to deal with these in this 
approach (see discussion in Ref. [28]). 

We now expand the resummed cross section, Eq. (fT4l) . in IPI kinematics through 
NNNLO. We provide results here for the case where Ts are trivial (1x1) color matri- 
ces. Explicit expressions for the more general case are found through NNLO in [24] and 
through NNNLO in [25]. 

At NLO, we find the expression for the soft-gluon corrections 

^(1) ^ ^B^Ml{c,Vi{s4)+C2Vo{s^)+Ci6{s4)} (17) 
vr 

where is the leading-order (LO) term, the LL coefficient is 

cs = Y,2a-Y,Cj, (18) 

i j 

with Cq = Cf and Cg = Ca, and the NLL coefficient C2 is defined by 02 = 02+^2, with 
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denoting the terms involving logarithms of the factorization scale, and 



T2 = 2Rer^s^-J2 

-E 



Ci + 2Ci In 



+ Ci\n 



M2 



5f + C,+C,ln^ 



(20) 



denoting the scale-independent terms. Again, ti denotes t or u, the sums over i run over 
incoming partons, and the sums over j run over any massless partons in the final state at 
LO. 

We write the NLO S{s4) terms as ci = + Ti, where 



Q. In 



denotes the terms involving logarithms of the factorization and renormalization scales. Here 
7^^^ = 3Ci?/4 and jg^^ = /3o/4, and Ti denotes virtual terms that cannot be derived from 
the resummation formaUsm but can be determined by matching to a full NLO calculation 
for any specified process. 

At NNLO, the soft-gluon corrections are 



(2) 



7r2 



\ cl P3(S4) + 



3 /?o , /^o 



+ 



C3Ci+C2-C2C3--r2 + -C3ln(^) +C3--^-5] 



+ 



(22) 



where we show explicitly results through next-to-next-to-leading logarithms (NNLL). For 
a complete expression see [24,25]. 

At NNNLO, the soft-gluon corrections are 



a 



a 



+ 



1 



I 4 2^5(54) + 



8 C3 C2 - — /3o C3 + — C3 /?0 2^ Cj 



Pi Po 



P4(S4) 



C3 C2 + 7^ C3 Cl - C2 C3 + — C3 - y C3C2 - YC3T2 + — C3 In I 



12 



6 



48 



^3(54) } + 



(23) 



where again we show explicitly results through NNLL. The complete expression is given 
in [25]. 
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In PIM kinematics we simply replace S4 by 1 — z, set s = M^, and drop the terms with 
ln(— tj/M^) in the above formulas. 

The NNNLO master equation, (l23l) . gives the structure of the soft corrections and can 
provide the full soft corrections explicitly if all the necessary two-loop and three-loop quan- 
tities are known. For processes with non-trivial color structure we are currently limited 
to NLL or NNLL accuracy. For processes with trivial color structure, such as 66 — > 
[27,29], all soft-gluon corrections have been determined through NNNLO. Below, the term 
means that the soft-gluon contributions through N^'^LL accuracy to the 
n-th order QCD corrections have been included. 

3 Applications of large-x resummations 

In this section we provide some calculations that are applications of the large-x resumma- 
tion formalism to processes of interest at the Tevatron and the LHC. We present results for 
top-antitop pair production, single top quark production, 1^-boson production, and Higgs 
boson production via 66 — > H. 

3.1 tt production 

The top quark, the heaviest known elementary particle, was discovered in pp collisions at 
Run I of the Tevatron in 1995 [30,31]. More recent measurements at Run II have increased 
the accuracy of the top mass and cross section measurements (for a review see [32, 33]) 
and thus require accurate theoretical calculations of top production cross sections and dif- 
ferential distributions. The main partonic channels in ti production are qq —>■ tt, which is 
dominant at the Tevatron, and gg ti, which will be dominant at the LHC. 

The latest calculation for top-antitop pair hadroproduction includes NNLO soft-gluon 
corrections to the double differential cross section [34]. Near threshold the soft-gluon cor- 
rections dominate the cross section at the Tevatron and contribute sizable enhancements. 
The form of the corrections and their numerical values depend crucially on the kinematics 
chosen to describe the process. The NNLO soft corrections were calculated fully to NNLL 
in both IPI and PIM kinematics [28, 35]. In addition a good approximation for the next- 
to-next-to-next-to-leading logarithms (NNNLL) was provided in [34]. The best theoretical 
result for the cross section is the average of the NNLO-NNNLL cross sections in the two 
different kinematics [34]. 

In Fig.[TJ we present the NLO and approximate NNLO-NNNLL tt cross sections at the 
Tevatron with = 1.96 TeV using the MRST2002 [36] parton densities. On the left we 
plot the cross sections as functions of mt, the top quark mass, for // = nit, where denotes 
the factorization and renormalization scales which we have set equal to each other. On the 
right we plot the cross sections as functions of fi/rrit with rrit = 175 GeV. The results are 
given in both IPI and PIM kinematics together with their average. The NLO cross sec- 
tion depends less on fi than the LO cross section, as expected. The NNLO-NNNLL cross 
sections exhibit even less dependence on fi, approaching the scale independence of a true 
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_ _ in — - 1/2 

pp->tt Tevatron S =1.96 TeV il=m^ PP > tt Tevatron S "=1.96 TeV m=175 GeV 




nij (GeV) n / 



Figure 1: The tt total cross sections in pp colhsions at the Tevatron with = 1.96 TeV 
as functions of rrit (left) and fi/mt (right). 

physical cross section. They change by less than 3% in the range mt/2 < ^ < 2mt. For a 
top mass of 175 GeV the average of the NNLO-NNNLL IPI and PIM results is 6.77 ± 0.42 
pb, where the uncertainty indicated is from the kinematics. Including all sources of uncer- 
tainty (kinematics, scale variation, and uncertainty from the parton distribution functions) 
we may write the cross section as 6.8 ± 0.6 pb. This theoretical result is in agreement with 
the latest experimental result for the cross section at the Tevatron [37,38]. Finally, we note 
that NNNLO soft-gluon contributions in the qq channel were presented in [25]. These 
NNNLO-NNLL corrections further stabilize the scale dependence of the cross section at 
the Tevatron. 

3.2 Single top quark production 

Single top quark production provides a way to directly measure electroweak properties 
of the top quark, such as the Vtt CKM matrix element. It also allows a deeper study of 
electroweak theory since the top quark mass is of the same order of magnitude as the elec- 
troweak symmetry breaking scale, and may be useful in the discovery of new physics. 
Therefore it is crucial to have accurate theoretical predictions for the cross section. 

The cross section for single top quark production is less than the tt cross section and 
the backgrounds to the production processes make the extraction of the single top signal 
challenging. Intensive searches for single top quark events at the Tevatron have recently 
produced evidence of such events [39,40]. The LHC has good potential for observation and 
further analysis of single top events. 

Single top quarks can be produced through three distinct partonic processes. One is 
t-channel production, qb —>■ q't and qb —>■ q't, via the exchange of a space-like W boson, a 
second is s-channel production, qcf — > bt, via the exchange of a time-like W boson, and a 
third is associated tW production, bg tW~. 

The threshold corrections to single top production have been calculated for both the 
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Figure 2: The cross section (left) and K factors (right) for single top quark production at 
the Tevatron in the s channel. Here n = fip = fiR = nit- 



Tevatron and LHC colliders through NNNLO [41-43]. At the Tevatron the t-channel pro- 
cess is numerically dominant, but the higher-order corrections are relatively small. The 
s-channel is smaller, but receives large corrections and it was shown that the threshold soft- 
gluon corrections dominate the cross section. Associated tW production is quite minor, 
although it also has large K factors, defined as the ratios of the higher-order cross sections 
to the LO cross section. At the LHC the t channel is again dominant, but the second largest 
channel is tW production; the s channel is numerically the smallest. Below we provide 
some numerical results for all three channels at both the Tevatron and the LHC colliders 
using the MRST2004 parton densities [44]. We add the soft-gluon corrections through 
NNNLO to the complete NLO cross section [45,46]. 

We begin with single top production at the Tevatron [41] with = 1.96 TeV. For 
t-channel production, the NNNLO-NLL cross section is cj*"'^'^^'^°'^^(mt = 175 GeV) = 
l.OS^OQ^ ± 0.06 pb, where the first uncertainty is from variation of the factorization 
and renormalization scales, fip and /ijj, between mt/2 and 2mt, and the second is due 
to the parton distribution functions. For the s channel, the corresponding cross section 
is cr'^-^hannoi^^^ ^ 175 GeV) = 0.49 ± 0.02 ± 0.01 pb. Finally, in the tW channel 
a*^{mt = 175 GeV) = 0.13 ± 0.02 ± 0.02 pb. We note that the cross sections for antitop 
production at the Tevatron are identical to those for single top production in each channel. 

In Fig. |2]we plot the cross section and the K factors for single top quark production at 
the Tevatron in the s channel setting both the factorization and renormalization scales to ^ = 
mt- We plot the LO cross section and the approximate NLO, NNLO, and NNNLO cross 
sections at NLL accuracy. The K factors are quite large, thus showing that the corrections 
provide a big enhancement to the cross section. 

We continue with single top production at the LHC [42] with \/S = 14 TeV. For 
the t channel the threshold corrections are not a good approximation of the complete 
corrections. The NLO cross section for top production al^p^^^'^^\mt = 175 GeV) = 
146 lb 4 lb 3 pb. For antitop production the corresponding result is o'antitop'^^^ ("^* ~ 
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b<;^->tW LHC S =14TeV 



bg-^>tW" LHC S =14TeV n=m 



LO 

■ NLO approx 

■ NNLO approx 
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NLO approx /LO 
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NNNLO approx / LO 




m (GeV) 
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Figure 3: The cross section (left) and K factors (right) for associated tW production at the 
LHC. Here ^ = = fJ-R = "T-t- 



175 GeV) = 89 ± 3 ± 2 pb. For the s channel, the soft-gluon corrections are rela- 
tively large and the soft-gluon approximation is good. The NNNLO-NLL cross section 
is at~p^^''''''\mt = 175 GeV) = 7.23t[3:^^ ± 0.13 pb for single top production and 
Knttor^^i^t = 175 GeV) = 4.03l[5i2 ± 0-10 Pb for single antitop production. Finally, 
for tW production the cross section is cr*^(mt = 175 GeV) = 41.1 it 4.1 it 1.0 pb, which 
is identical to that for associated antitop production. In Fig. [3]we plot the cross section and 
K factors for associated tW production at the LHC setting the scales to ^ = As seen 
from the plots, the soft-gluon corrections are large for this process. 



3.3 l^-boson production at large transverse momentum 

VF-boson production in hadron colliders can be used in testing the Standard Model and 
in estimating backgrounds to Higgs production and new physics. Precise calculations for 
W production at large transverse momentum, Qt, are needed to identify signals of new 
physics which may be expected to enhance the Qt distribution at high Qt- 

Analytical NLO calculations of the cross section for W production at large transverse 
momentum were presented in Refs. [47,48], where numerical results were also presented 
for the Fermilab Tevatron. Numerical NLO results for W production at the LHC were 
more recently presented in [49]. The NLO corrections enhance the Qt distribution of the 
W boson and they reduce the factorization and renormalization scale dependence of the 
cross section. 

A recent theoretical study [50] included soft-gluon corrections through NNLO, which 
provide additional enhancements and a further reduction of the scale dependence. The 
complete NNLL terms were calculated and an approximation for the NNNLL terms was 
derived at NNLO. Numerical results with these soft corrections were calculated for W 
production at the Tevatron [50] and the LHC [49]. 

Here we discuss W production at large transverse momentum at the LHC with = 14 
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TeV using the MRST2002 parton densities [36]. The LO partonic processes for the produc- 
tion of nW boson and a parton are qg — > Wq and qq Wg. The electroweak coupling 
a(M|) is evaluated at the mass of the Z boson, and standard values [51] ai^e used for the 
various electroweak parameters. In the numerical results we present the sum of cross sec- 
tions for and production. The W bosons at the LHC will be detected primarily 
through their leptonic decay products e.g., W~ — > Iv^, therefore the cross sections pre- 
sented here should be multiplied by the appropriate branching ratios. 

In Fig. |4](left plot) we plot the transverse momentum distribution, da/dQ^, at high Qt 
for W production at the LHC. We set fip = fiR = Qt and denote this common scale by 
fi. We plot LO, NLO, and NNLO-NNNLL results using the corresponding parton densities. 
As seen from the plot, the NLO corrections provide a significant enhancement of the LO Qt 
distribution. The NNLO-NNNLL connections provide a rather small further enhancement of 
the Qt distribution. However, the NNLO-NNNLL corrections can be much bigger for other 
choices of factorization and renormalization scales. The NLO corrections increase the LO 
result by about 30% to 50% in the Qt range shown. In contrast, the NNLO-NNNLL/NLO 
ratio for this scale is rather small. Part of the reason for this is that the NNLO parton 
distribution functions are significantly smaller than the NLO pdf. 

On the plot on the right in Fig. |4]we show the scale dependence of da/dQj^ for Qt = 
150 GeV versus ^i/Qt over two orders of magnitude. It is interesting to note that the scale 
dependence of the cross section is not reduced when the NLO corrections are included. 
This is due to the fact that the cross section is dominated by the process qg —>■ Wq. The 
gluon density in the proton, at fixed x less than '^O.Ol, increases rapidly with scale. Thus, 
the fiji and fip dependencies cancel one another to a large extent. However, we have an 
improvement in the scale variation when the NNLO-NNNLL corrections are added. The 
NNLO-NNNLL result displays very little scale dependence. 
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bb->H at Tevatron S "=1.96 TeV n=m 
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Figure 5: The K factors for 66 — > at the Tevatron (left) and the LHC (right). Here 



3.4 Higgs boson production via bb H 

The search for the Higgs boson [52] is one of the most important goals at the Tevatron and 
the LHC colliders [53]. The main Standard Model production channel at these colliders is 
gg H. However, the channel bb H can be competitive in the Minimal Supersym- 
metric Standard Model at high tan /3, with tan (3 the ratio of the vacuum expectation values 
for the two Higgs doublets. The complete NNLO QCD corrections for this process were 
calculated in [54]. 

Complete expressions for the soft-gluon corrections at NNNLO were presented in 
[27, 29]. However, it is known at NNLO that the soft corrections alone are not a good 
approximation of the full corrections [27,54]. Purely coUinear terms [27,55,56] have to 
be included to provide an accurate calculation. An approximation for the coUinear terms 
through NNLL accuracy at NNNLO was provided in [27]. 

We now present numerical results for 66 ^ at the Tevatron and the LHC [27] using 
the MRST2006 parton densities [57]. Figure |5] shows the K factors for Higgs production 
via bb ^ H at the Tevatron (left) and the LHC (right), with fj, = mu- The complete 
NLO corrections increase the LO result by around 60% at both the Tevatron and the LHC. 
Inclusion of the complete NNLO corrections futher increases the cross section: the NNLO 
K factor is around 1.9 at the Tevatron and 1.8 at the LHC. By including at NNNLO the sum 
of the complete soft-gluon corrections and the coUinear approximate NNLL corrections 
(S+NNLLapp), we find further enhancement. From the study of the contributions of the soft 
and coUinear terms at NLO and NNLO at both the Tevatron and the LHC we expect that 
the NNNLO S+NNLCapp curve provides a good approximation of the complete NNNLO 
cross section. The NNNLO S+NNLCapp K factor is between 2.06 and 2.01 at the Tevatron 
and between 1.95 and 1.87 at the LHC for Higgs masses ranging between 110 and 180 GeV, 
which is a significant addition to the NNLO result. 
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4 Loop calculations in the eikonal approximation 

The soft-gluon resummation formalism, and in particular the calculation of the soft anoma- 
lous dimension matrices, employs the use of the eikonal approximation in loop diagrams. 
The eikonal approximation is valid for descibing the emission of soft gluons from partons in 
the hard scattering. The approximation leads to a simplified form of the Feynman rules by 
removing the Dirac matrices from the calculation. When the gluon momentum goes to zero, 
the Feynman rules for the quark propagator and quark-gluon vertex in Figure [6] simplify as 
follows: 

u{p) i-igsTp) 1^^—7X2 ^ ^yP> 9sTp Tt, 7—- = ^yP) 9sTp — —— 

[p + fcj^ — + i€ 2p ■ K + le V ■ k + le 

(24) 

with V a dimensionless vector, p c(. v, and Tp the generators of SU(3) in the fundamental 
representation. 



p + k 1 


Z P 







Figure 6: Eikonal approximation. 

The ultraviolet poles in loop diagrams involving eikonal lines are particularly important 
as they play a direct role in the renormalization group evolution equations that are used in 
threshold resummations [4,6,7] (see Eq. (fT3])). 

Below we give examples of a one-loop and a two-loop calculation for diagrams involv- 
ing eikonal lines representing massive quarks. For the calculation we use the Feynman 
gauge, and we use dimensional regularization with n = 4 — e dimensions. 

4.1 One-loop calculation 

In this subsection we calculate the integral In for the one-loop diagram in Fig. ^ with 
eikonal lines representing outgoing massive quarks. This one-loop integral is given by 

^' J (27r)" A;2 „ . . ^ . ^) ■ ^^^> 
Using Feynman parameterization, this integral can be rewritten as 

hi = -Ws 7^ \ dx\ dy — — -— ^ -3 . (26) 

(ivr)" Jo Jo J [xk^ + yvi ■ k + {I - X - ii)vj ■ k] 
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Figure 7: One-loop eikonal diagram with outgoing massive quarks. 



After several manipulations, Eq. (1261 ) becomes 

hi = ^ (-1)-W2 25^/2 ^^/2 r (^1 + (l+/32)^'^^^-l+e(;L_^)-l-e 

^ {/; - K<^ - ^) - "T ' - f r -^^^ife^ - ° ) 

(27) 

where here /5 = \/l — Arri-^ /s, with m the quark mass, and we have used the relations 
Vi ■ Vj = (1 + /32)/2 and v^ = v] = (1 - P^)/2. 

The integral over x in Eq. (l27l) contains both ultraviolet (UV) and infrared (IR) singu- 
larities. We isolate the UV singularities and find that 



rl 1 

/ dxx-^+'{i-x)~^-' = - + m. 

Jo e 



(28) 

After calculating the integrals over z in Eq. (1271 ). we find that the UV poles and constant 
terms of In are 

/[[^ = ±ljlJ_ln ^ +_(41n2 + ln7r-7£;-z7r)ln' 



vr 2/3 ie Vl + Z?/ 2' ' ' \1 + P . 

.l>„^a.«-i>n^a-«-iL.(l±i?).iL.(l^)}.,„ 

Complete one-loop calculations for heavy quark production in axial gauge were pre- 
sented in Ref. [4]. 

4.2 Two-loop calculation 

In this subsection we calculate the two-loop integral I21 for the quark-loop diagram in 
Fig. ([Hi given by 



hi = i-^hfgt J 



d^k d^l V>^ i-V^) {-i)g^,{-i)g 

(27r)" (27r)" Vi ■ k {-vj • k) k"^ k"^ 
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Figure 8: Two-loop eikonal diagram, involving a quark loop, with outgoing massive quarks. 



xTr 



After a few manipulations involving the trace we can write this integral as 

■,4 



^21 



-Ant 



9 s 



where 



f (27r)2« 
I21 = ■ Vj 

= -Vi ■ Vj 



^21 + ^21 + ^21 + ^21 + 



kvj-kk'^ J {I- k) 



dTk 



/ • k 



kvj-kk^ J P{1 - k) 



'2/ 



Vi-kVj-kk^j P{l-k)^ 
d-k 
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2 * k J\i^ 
Vj ■ k A;4 



P{1 - kf 
P{1 - kf ■ 



We begin with the evaluation of I^i- Since 



{I - kf 



we find If, = 0. 



Next we evaluate I2;. Using Feynman parameterization, we find 



jd-l 



P{1 - k) 



(30) 



(31) 

(32) 
(33) 
(34) 
(35) 
(36) 

(37) 



(38) 
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therefore 



'2« 



-iir 



(5-e)/2 2-2+. ^..^.p(e 



r(i 



(39) 



The k integral in the above expression is 

d''k 



Vi-kvj-k (A;2)i+^/2 



ITT 



2-e/2 22+2e (^_iyl-e 



r(i + 



X / dxx 



-l+2e 



(1-x) 



-l-2e 



1-/32 



n -l-e 



.(40) 



The integral over x in Eq. (|40l l contains both UV and IR singularities. We isolate the UV 
singularities and find that 



Jo 



dx X 



-l+2e 



{i-xy 



-l-2e 



1 

Ye 



+ IR. 



(41) 



The integral over y is given in terms of hypergeometric functions 
2F1 ^— e, 1 + e, 1 — e, which can be expanded in powers of e. After some cal- 



culation, we find 

d'^k 



Vi-kvj-k {k^y+^/^ 



2m^ 

~p7 



In 



+ 



2m^ 



Li2 



ln2(l + /?) - ln2(l - /?) + -(6In 2 - In^ - 7ij) In 



1 + 
1-/? 
1 + /? 



Li2 



1-/3. 
(42) 



Assembling everything together we find the result for the UV poles of 



jh UV 



IT 



.(1+^ 

(3 



1-/? 
1 + /? 



1 

+ - 

e 



Lio 



1 + 



Lio 



+ ln2(l ln2(l - /?) + (1 + 31n2 - Invr - 7s) In 



1-/3 
1 + /3 



.(43) 



We continue with the evaluation of Igi- Now 



/2(Z - kf 



m2-^/2r 

,2-e/2 



r 1 



r(3-f) 

2j r(4-e) 



ITT 



+ 



and after a few manipulations we find 



-g^'T -1 + 



^A(r(2-f))^ 



r(4 - 6) 



{k 



2\l-€/2 



(44) 



-'2Z 



Vi-kvj-k (A;2) 1+^/2 



(45) 
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The integral over k was evaluated before for Eq. (l40l ). We thus find that the UV poles 



of Ifo 



TCUV 

hi 



vr 



3/3 



■In 



1 + /? 



+ 



1 



-Li2 



1 + /3 



+ Li2 



ln"(l + /3) + ln"(l - /3) + -- - 31n2 + Invr + In 



1-/? 
1-/? 
1 + /3 



(46) 



Finally, we calculate and /|;. We use Eq. (|38] ) for the / integral and then find that 
the remaining integral over k vanishes, so = I|j = 0. 

Adding all the terms in Eq. (|3TI ). the final result for the UV poles of I21 is 



tUV 

J- 21 



-n 



a2(l + /32) r 1 /l-/3\ , 1 



^2 6/3 



In 



1 + /3 



+ 



Li2 



1 + /3 



+ ln^(l + /?) - ln^(l -/?)+(- + 51n2 + In vr - 7£; ) In 



- Li2 

1-/3 
1 + /3 



1-/3 



(47) 



More results for two-loop integrals with massive quarks wiU appear in [23] . 



5 Small-x resummations 



In the last forty years there has been a large effort trying to understand what are the correct 
effective degrees of freedom underlying the strong interaction at high energies. In scat- 
tering processes where the center-of-mass energy is much larger than any other scales the 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) approach [10-14] emerges as the correct approach 
to describe the scattering. This framework relies upon t-channel "Reggeized" gluons in- 
teracting with each other via standard gluons in the s-channel and a gauge invariant three 
particle vertex. This simple structure is a consequence of using multi-Regge kinematics 
where gluon cascades are ordered in longitudinal components but with a random walk in 
transverse momenta. Although this simple iterative and linear structure must be modified 
at higher energies in order to introduce unitarization and non-linear corrections, there is a 
window at present and future colliders where the BFKL predictions hold. 

In the leading logarithmic approximation (LLA) we resum terms of the form (a^ In s)". 
Diagrams contributing to the running of the strong coupling do not appear and the coupling 
is a constant parameter. The factor needed to scale the energy in the logarithms is also free 
and the predictability of the LL approximation is limited. In the next-to-leading logarithmic 
approximation (NLLA) diagrams with an extra power in the coupling without introducing 
an extra logarithm in energy are considered. The coupling is allowed to run and the energy 
scale is determined. 

In this contribution we discuss three aspects of the BFKL resummation program. In 
subsection 15.11 we review the relevant equations to describe final states at small values of 
Bjorken x in Deep Inelastic Scattering (DIS). We introduce the concept of color coherence 
and the CCFM equation. We show the differences and similarities between the BFKL 
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approach and the introduction of angular ordering in the case of jet rates. In subsection 15 .21 
we analyse in detail how to extend the region of applicability of the multi-Regge kinematics, 
the basic ingredient in the BFKL approach, to regions with coUinear emissions. We will 
find an interesting structure in the higher-order corrections that can be resummed into a 
Bessel function of the first kind, which accounts for the double logarithms in tranverse 
scales. In subsection [53] we briefly explain the SL{2, C) invariance associated to the BFKL 
Hamiltonian and how it shows up in the physics of multijet events, in particular in the 
production of Mueller-Navelet jets at a hadron collider. 



5.1 QCD coherence and small-x final states 

In Quantum Electrodynamics coherence effects are responsible for the suppression of soft 
bremsstrahlung from electron-positron pairs. In QCD processes such as g qq any soft 
gluon emitted with an angle from one of the fermionic lines larger than the angle of emission 
in the qq pair will probe the total color charge of the pair. This charge is the same as the one 
from the parent gluon and the radiation takes place as if the soft gluon was emitted from it. 
This color coherence leads to the angular ordering of sequential gluon emissions. 

In DIS, let us say that the (i — l)th emitted gluon from the proton has energy Ei^i. A 
gluon radiated from it with a fraction (1 — Zj) of its energy and a transverse momentum 
will have an opening angle 



« (48) 

(1 - Zi)Ei.i 



with 



p - (49) 

Color coherence leads to angular ordering with increasing opening angles towards the hard 
scale (the photon). Therefore, we have O-i^i > 9i, or 

Q.+i >_£l^, (50) 



1 - Zi+l 1- Zi 

which in the limit Zi, Zj+i <^ 1 reduces to 

qi+i > Ziqi. (51) 

In Ref. [58-61] the BFKL equation for the unintegrated structure function was obtained in 
a form suitable for the study of exclusive observables: 

Uk) = 1^ t -z-^r{z, k)e{q - fi)Uq + k), (52) 

J TTq^^ Jo z 

where // is a collinear cutoff, q is the transverse momentum of the emitted gluon, and the 
gluon Regge factor is 



AR{zi,ki) = exp 



- 1 1 , 
-as m — m — 

Zi 11^ 



(53) 
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with ki = \ki\, and as = asNc/ir. Under iteration, this expression generates real gluon 
emissions with all the virtual corrections summed to all orders. Since is an inclusive 
structure function, it includes the sum over all final states and the /x-dependence cancels 
between the real and virtual contributions. 

The structure function is defined by integrating over all < qf < Q^, i.e. 

Fo,(Q,Ai) = e(Q-/x) + 5] / \{-Sldzi'^ztl^R{zi,h), (54) 

with i real gluon emissions in each iteration of the kernel. The contributions from a fixed 
number r of emitted gluons is 



/■I °° 
Fo^iff) = / dx x^Foix, Q) = 1 + ^ fI\q). 

In Ref. [59] the perturbative expansion for the F^J {Q, /j,) 

oo —n 

Fi:!{Q,f^) = T.ci[\n;T)'^, 

n=r 

was obtained with T = ln{Q/iJ,). Then we have 



(55) 



(56) 



FoUQ)^T.^SiQ) 



i=Q 



Q_ 



2\ ^ 



(57) 



where 7 is the BFKL anomalous dimension. It was pointed out that coherence effects sig- 
nificantly modify the individual Fq'J (Q) whilst preserving the sum Fouj{Q), and care must 
be taken to account properly for coherence in the calculation of associated distributions. 
Modifying the BFKL formalism to account for coherence [58-61], Fo^^{Q, fi) becomes 

00 „q2 t i2 

FUQ,fi) = e{Q-fi) + Y^ / l[-^dz,^zfA{zi,qi,ki)@{qi-Zi-iqi-i), (58) 



r=l 



where AR{zi,ki) is substituted by 

A{zi,qi,ki) = exp 



=1 



- , 1 , 
-as in — m — ^ 



; ki > Qi, 



and for the first emission we take qozo = fi. The expansion of F^^^ (Q) is now 

00 n -n 



(59) 



n=r m=l 



(60) 



A coUinear cutoff is needed only on the emission of the first gluon because subsequent 
coUinear emissions are regulated by the angular ordering constraint. 
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The rates for the emission of a fixed number of resolved gluons, with a transverse mo- 
mentum larger than a resolution scale jiR, together with any number of unresolved ones, 
were calculated in Ref. [62] in the LLA, to third order in as- fiR is constrained by the 
coUinear cutoff and the hard scale, ii <^ ijlr <^ Q. For the n-jet rate all the graphs with 
n resolved gluons and any number of unresolvable ones were considered. Expanding the 
Regge factors to 0(05) we find that the jet rates both in the multi-Regge (BFKL) approach 
and in the coherent (CCFM) approach are the same: 



Ojet 

Ijet 

2 jet 

3 jet 



{2as)^ ^ {2asf 
{2as)^ , {2asf 



52 



+ 



{2asf 



T + 



OJ 

{2as? ^ 
0? 

[2asf 



OJ^ 



TS - -T^ 



+ 



+ 
7 



53 
6 

(205) 



to-' 



3 



2 2 



(61) 

,(62) 
(63) 
(64) 



with T = \n{Q/fiji) and S = ^). When coherence is introduced the singularities at 

(J ~ are stronger than in the BFKL approach but the extra logarithms cancel in the sum of 
all the graphs needed for the jet rates. The net effect is that the final results are the same as 
those obtained without coherence [62]. This is true to all orders in the couphng [63] since 
a generating function for the jet multipUcity distribution was obtained in [64]. Within the 
multi-Regge kinematics the r-jet rate reads 



1 a" 



where the jet-rate generating function R^j is given by 



(65) 



u=0 



2a 



1 + (1_^)^T 

OJ 



(66) 



The same generating function is obtained when coherence is considered. The mean number 
of jets and the mean square fluctuation in this number are 



(n) 



du 



RUu,T) 



u=l 



LO 2 \ UJ 



(67) 



,2. / ^2 2a„ 3 /2a. \2 2 /2a. \3 

00 2 \ OJ J 6 \ OJ J 

In general, the pth central moment of the jet multiplicity distribution is a polynomial in 
OLgT joj of degree 2p — 1, indicating that the distribution becomes relatively narrow in the 
hmit of very small x and large Q/ ijlr [64]. 
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In Ref. [65, 66] the subject was developed even further and all subleading logarithms 
of Q'^/^J-ji were included to calculate the jet multiplicity in Higgs production at the LHC. 
In Ref. [65] they extended the results from a [a5ln(l/a;) ln{Q'^ / fij^)]"- resummation to a 
[0:5 ln(l/x)]"[ln((3^///^)]'" one with < m < n, proving that the quadratic and cubic 
forms of the mean and the variance remain vaUd. It has also been shown that for any suf- 
ficiently inclusive observables the CCFM formalism leads to the same results as the BFKL 
equation [67]. The key idea to understand this result comes if we try to obtain the results of 
Ref. [59] from those in the previous section in the limit i^r 0. To get the right solution 
we should consider subleading terms ocg ^t^^{Q/i^r) which must be resummed when taking 
the limit fiR — > to obtain a continuous transition from the case where BFKL and CCFM 
results are equivalent, to that of them being different. If we also consider the effects of 
introducing the 2; ^ 1 divergent part of the splitting function in the CCFM approach we 
will see that this leads to all BFKL and CCFM final-state properties being identical in the 
[qs ln(l/x) ln(Q^//i|,)]" approximation [67]. Recent reviews devoted to the implementa- 
tion of CCFM in Monte Carlo event generators can be found in, e.g., [68-71]. An approach 
which has the potential to apply BFKL in the NLLA to DIS phenomenology is that in 
Ref. [72-75]. In the NLLA approximation it is important to carefully take into account kt 
factorization [76]. 



5.2 Improving the coUinear region of multi-Regge kinematics 



In this section we revisit the approach of Ref. [77] where the multi-Regge kinematics was 
extended to include coUinear contributions present to all orders in the BFKL formalism. In 
Ref. [78] the structure in transverse momentum space of the double logarithms resummed 
was exphcitly extracted. A new renormalization group (RG)-improved kernel was obtained 
which does not mix transverse with longitudinal momentum components. 

In the MS renormalisation scheme, the BFKL kernel in the NLLA acting on a smooth 
function [15, 16] is 



I 



/3o 



In 



2|2 
92 I 



ANc \ max {ql,q^) fi^ 



/(<,!)- 



' {Q1+Q2) 



f(<^^))-^(r{qu^)+ln■,^2 





(69) 



where (3o = {llNc - 2n/) /3, S = (4 - tt^ + 5/3o /iVc) /12, and T{qf, ql) can be found 
in Ref. [15]. The coUinear structure can be obtained acting on the eigenf unctions in the 
LLA, i.e. 



(fq2)C{qi,q2) 



' as {ql) 
as {ql). 




as {ql) Xo (7) + alxi (7) • (70) 
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Here we have 

Xo(7) = 2V(1) - (7) - V' (1 - 7) , (71) 

XI (7) = 5x0 (7) + ^ « (7) +V'" (1-7)) -^(0(7) + 0(1 -7)) (72) 

n'cosjnj) / I I (2 + 37(1- 7)) \ , 3 _ ^ 2 . . 

4sm2(7r7)(l - 27) V I J (3 - 27) (1 + 27);' ^ 2''' 8iVe^° ' 

with V (7) = r' (7) /r (7) and 

.^(7)+<^(l-7) = 



m=0 



The pole structure around 7 = 0, 1 is 



Xo(7) - - + {7^1-7}, (74) 

7 



XI (7) - - + ^- A + {7^1-7} (75) 

7 7^ 27"^ 



with 



5_Po_ _ 13 n/ _ 55 _ 1^ _ _ 11 

12 AT^ 36A^3 36' g jy^ 6iV3 ^ ^ 

The cubic poles stem from t/j" and compensate for the equivalent terms appearing when 
the Regge-like energy scale sq = qiq2 is shifted to the DIS choice sq = Qi,2- Higher 
order terms beyond the NLLA, not compatible with RG evolution, are also generated by 
this change of scale. The NLLA truncation of the perturbative expansion is then the reason 
why the gluon Green's function develops oscillations, where the Green's function can have 
negative values, in the qf / ratio. 

It is possible to remove the most dominant poles in 7-space incompatible with RG 
evolution by simply shifting the cu-pole present in the BFKL scale invariant eigenfunction. 
Here we focus on the scheme proposed in Ref. [77]: 



a; = a, I 1+ I a+^j Q,j (^2^(1) - ^ (^7 + | " b a, j - ^(^l-7 + |- ba, 



+ a' [xi (7) + (^Xo (7) - b) (^'(7) + V''(l - 7)) - + y) ^0(7)) • (77) 
We can approximately solve this equation considering the w-shift in the form 

(1 + Aa,) (^2V'(1)-V(7 + | + Ba,) - V (^1 - 7 + | + Ba.)) , (78) 



uj = a. 
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which can be written as 



as (1 + Aa,) ^ 



+ 



m=0 



(79) 

The solution to this shift can be obtained by adding all the approximated solutions at the 
different poles plus a term related to the virtual contributions, i.e. 



OJ 




(1 + 2m + 2Bas) + |7 + m + Ba^ 
2as (1 + Aa,) 



1 + 



2as (1 + kas 



+ 1-7 + m + Bas 1+ ^ 

'V (1-7 + m + Ba,)^ 

At the 7 = 0,1 poles this expansion generates the NLLA terms: 



(7 + m + B as 

2as (1 + Kas 
m + 1 



(80) 



a. 



UJ 



7 



+ {7^1-7} 



'A B 1 

7 72 27*^ 

To match the original kernel at NLLA we set A = a and B = — b from Eq. ( [76l ). 

We now include the full NLLA scale invariant kernel without double counting terms: 



(81) 



(82) 



+ 




(-l)"(2n)! (a, + aag) 

\2n+l 



-Q2"n!(n + 1)! + m 
a b 



b a. 



a. 



+ 



1 



+ {7^1-7} 



7 + m " \'j + m ' (7 + m)^ 2(7 + m)'^. 

This result reproduces the w-shift very closely, see Fig. |9] The imaginary part of 7 at the 
maximum of the NLLA scale invariant eigenvalue (middle plot of Fig. |9l) is not zero and 
results in oscillations in the qf / variable. These are eliminated when the RG-improved 
kernel is used, as it also happens for the "all-poles" kernel. 

It is very important to note that in Eq. (l82l ) the tij-space is decoupled from the 7- 
representation. In Ref. [78] an expression for the coUinearly improved BFKL kernel which 
does not mix longitudinal with transverse degrees of freedom was found. The only mod- 
ification needed in the full NLLA kernel to introduce the "all-poles" resummation is to 
remove the term 



a: 



1 



in the real emission kernel, K,r (^q,kj, and replace it with 
1 



(83) 



(q- kf 




k — q 



\ 



2 (tts + aa'i] 



Ji W2(q, + aa2)ln^ 
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Figure 9: 7-representation of the LLA and NLLA kernels. The RG-improved kernel by a 
cj-shift is included together with the new "all-poles" approximation. 



-a,-aa^+ba^L_^ln(^|^j|, (84) 

with Ji the Bessel function of the first kind. When the difference between the and k'^ 
scales is not very large then 




and its influence is minimal, not affecting the "Regge-Uke" region. When the ratio of 
transverse momenta becomes larger then 

(f . (g)) = [-1-^^ ' - {() - f ) 

compensating for the unphysical oscillations. This resummation of all-poles has been ap- 
plied to extend the region of applicability of BFKL calculations in the NLLA in the case of 
electroproduction of Ught vector mesons in Ref. [79]. 

5.3 Conformal signatures at the Large Hadron Collider: azimuthal angle 

We now proceed to review the work of Ref. [80] where azimuthal angle decorrelations in 
inclusive dijet cross sections were studied analytically to include the NLLA to the BFKL 
kernel, while keeping the jet vertices at leading order. It was shown how the angular decor- 
relation for jets with a wide relative separation in rapidity largely decreases when higher 
order effects are considered. 

Observables where BFKL effects should be dominant require a large enough center- 
of-mass energy, and two large and similar transverse scales. An example is the inclusive 
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hadroproduction of two jets with large and similar transverse momenta and a large relative 
separation in rapidity, Y, the so-called Mueller-Navelet jets, first proposed in Ref. [81]. A 
rise with Y in the partonic cross section was predicted in agreement with the LLA hard 
Pomeron intercept. At the hadronic level, Mueller-Navelet jets are produced in a region 
where the parton distribution falls very quickly, reducing this rise. Small x resummation 
effects are very relevant if we investigate the azimuthal angle decorrelation of the pair of 
jets. BFKL enhances soft real emission as Y increases, reducing the angular correlation. 
This was first investigated in the LLA in Ref. [82-84]. The rate of decorrelation in the LLA 
lies quite below the experimental data [85-88] at the Tevatron. This motivates the NLLA 
discussion of this subsection. 

We are interested in the cross section parton -i- parton jet -i- jet -i- soft emission, with 
the two jets having transverse momenta and q2 and with a relative rapidity separation Y. 
The differential partonic cross section is 

^ 7r^«3/(gi,g2,Y) ^^^^ 

It is useful to introduce a Mellin transform: 

fiqi,q2,Y) = /|^e"'^/-(9"'i'^2). (88) 
The solution to the BFKL equation in the LLA is 

in{ei-e2) 



Xo{\n\,iy) 



(89) 



with 



Xo{n,i^) = 2V'(l)-vQ+i'^ + ^) -V'Q-^^+^)- (90) 

The nonforward BFKL equation corresponds to a Schrodinger-like equation with a holo- 
morphically separable Hamiltonian where — i Y is the time variable. Both the holomorphic 
and antiholomorphic sectors are invariant under spin zero Mobius transformations with 
eigenf unctions carrying a conformal weight of the form 7 = ^ + iv + ^. In the principal 
series of the unitary representation, v is real and |n| the integer conformal spin [89]. In this 
way extracting information about n is equivalent to proving the conformal structure of high 
energy QCD. 

We now integrate over the phase space of the two emitted gluons together with some 
general jet vertices, i.e. 

In the jet vertices only leading-order terms are kept: 

<U^>?) = (^i'i'-pf)^ (92) 
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where pf corresponds to a resolution scale for the transverse momentum of the gluon jet. 
To extend this analysis it is needed to use the NLO jet vertices in Ref. [90, 91] where the 
definition of a jet is much more complex than Eq. ( |92l ). We then have 



f{qi,q2,Y). (93) 



In a transverse momenta operator representation: 



ir)' 



"2 inO 



(94) 



the action of the NLO kernel, calculated in Ref. [92], is 



K\u,n) = 

-2 PO 



as Xo {\n\ ,1') + alxi {\n\ ,v) 



d 



2X0 {\n\,y) ( +log^ ) + ( i—Xo 



d 



di 



W,n), (95) 



where xi, for a general conformal spin, reads 



Xi(n,7) = '5xo(?^,7) + ^C(3) - ^Xo(">7) 



+ 



1 



vr^ cos {-Kj) 



4sin^ (7r7) (1 - 27) 



o I ^ I 2 + 37(1-7) 

. ^ N^J (3-27)(l + 27) 



V N^J 2 (3 -27) (1 + 27)°"' 



(96) 



The function (/> can be found in Ref. [92]. 

The jet vertices on the basis in Eq. (l94l) are: 



/ 



A" 



^n,0 = Cl (z^) 5n,0, 



being replaced by pg- The cross section can then be rewritten as 



(97) 



, 1 1 

— {q \v,n) 

with the C2 iy) projection of ^f^^ on (n, i/| g) being the complex conjugate of ( [97] ) with p\ 



a 



Pl,2 



7r2a2 °° 



— „ J — CO 



X U + a. Y 



AN, 



XI i\n\ , i^) + TTT log (/X ) + -— lo; 



Cl (z^j 



2 9z^ \C2{iy)J 2du 



(98) 



+ ) XO (|n| ,1^) 



Resummations in QCD hard-scattering at large and small x 



29 




Figure 10: Evolution of the partonic cross section with the rapidity separation of the dijets. 



For the LO jet vertices the logarithmic derivative in Eq. (1981) is 

= log (phi) + 



dv \C2{y) 



(99) 



If (/) = ^1 — ^2 — TT, in the case of two equal resolution momenta, pf = P2 = P^' angular 
differential cross section can be expressed as 



with 



C„(Y) 



" ^ n.= — nn 



«.(p^)YL(|n|,.)+a,(p2)L{M,.)-J^^il(M^ 

du e V \ \^ " ) 



(100) 



(101) 



The coefficient governing the energy dependence of the cross section corresponds to n = 0: 

„2^ 



P 



2p2 



Co(Y). 



(102) 



We have chosen the resolution scale p = 30 GeV, = 4 and Aqcd = 0.1416 GeV. The 
71 = coefficient is directly related to the normalized cross section 



a(0) 



Co (0) ■ 



(103) 



The rise with Y of this observable is shown in Fig. [TO] Clearly the NLL intercept is 
very much reduced with respect to the LL case. The remaining coefficients with n > 1 
all decrease with Y. Because of this, the angular correlations also diminish as the rapidity 
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Figure 11: Dijet azimuthal angle decorrelation as a function of their separation in rapidity. 

interval between the jets gets larger. This point can be studied in detail using the mean 
values 



(cos {m(j))) 



Cm{Y) 



(104) 



Co(Y)- 

(cos {(/))) is calculated in Fig. [TT] The NLL effects decrease the azimuthal angle decor- 
relation. This is the case for the running of the coupling and also for the scale invariant 
terms. This is encouraging from the phenomenological point of view given that the data at 
the Tevatron typically have lower decorrelation than predicted by LLA BFKL or LLA with 
running coupling. The difference in the decorrelation between LLA and NLL A is driven by 
the n = conformal spin since the ratio 



(cos (0)) 



NLLA 



^NLLA (Y) Cj^LA 



(Y) 



(cos {(p)) 



LLA 



rNLLA 



(Y) C}^LA (Y) ' 



remains in the region 



1.2 > 



(^NLLA 



(Y) 



(Y) 



> 1. 



(105) 



(106) 



This is a consequence of the good convergence, in terms of asymptotic intercepts of the 
NLLA BFKL calculation, for conformal spins larger than zero. For completeness the 
m = 2,3 cases for (cos {mcj))) are shown in Fig. [12] These distributions test the struc- 
ture of the higher conformal spins. The methods of this subsection have been applied to 
phenomenology of dijets at the Tevatron and the LHC in [93,94], and to the production of 
forward jets in DIS at HERA in [95]. 



6 Conclusion 



The precision of perturbative QCD calculations will play a major role in the confidence of 
new physics discoveries, both at this generation of experiments, Tevatron and LHC, and 
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in future experiments. The most available avenue of improving the precision of QCD is 
through resummation of large contributions. We have presented results for the resumma- 
tion of large-x contributions and separately small-x contributions. In both cases, the large 
contributions arise from incomplete cancellations of virtual and real terms, and can be com- 
puted in the eikonal approximation. 

We have shown that the inclusion of soft-gluon corrections to top quark production 
cross sections is essential to stabilize the unphysical scale variations in the order-by-order 
calculations. This is necessary for any sort of precision calculation of the top mass and 
production channels. Additionally, we have shown the importance of resummation on W 
production at large transverse momentum, and on Higgs production. Discovery of the Higgs 
boson is the last remaining test of the Standard Model and precision measurements of its 
properties is essential to proceed forward with beyond the Standard Model theories. 

We have also presented a framework to include colUnear effects into the BFKL formal- 
ism. This stabilizes the oscillatory behavior that arises when one moves away from the strict 
kinematic regime of validity. It was shown how this inclusion improves the prediction of 
Mueller-Navalet jets, jets with a large rapidity separation but similar transverse scales. This 
is a process which will be observed at the LHC where the BFKL formalism should flourish; 
an important test of the complex behavior of QCD. A comparison between the predictions 
steming from a pure BFKL analysis and one including QCD coherence in multijet final 
states in DIS has been also discussed in detail. 
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